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ABSTRACT 


Capacitance and conductance measurements using two and 
four wire techniques were developed and statically and 
dynamically calibrated in this thesis. The voltage 
sensitivities range from 7.3 to 8.1 += 0.1 mV/cm for the ERA 
wire capacitance system static calibrations. This is ż 5.2% 
of the limiting theoretical value. The voltage sensitivities 
range from 0.3 to 0.4 + 0.1 V/cm for the wtour aa 
conductance system static calibrations. Dynamic calibrations 
were only completed for the conductance system. The dynamic 
calibration results were weakly frequency dependent with a 
170.15 decay in a limited, 2-4 Hz range. 

Wind power spectrum measurements were taken in the 
existing Upper Ocean Simulations Facility at the Nazis 
Postgraduate School. There was excellent agreement in the 
spectra with both techniques. 

Driven gravity wave frequency downshifting and wind 
energy dumping was observed in the combined gravity wave and 


wind-wave measurements. The power spectra peaked near two 


Hertz and decayed at 50 to 70 dB per decade, or as 073 to 
07? for both systems. 

Gravity wave phase speed and wavelength measurements 
were performed with the capacitance system. The results were 


approximately 40% higher than theory. 
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I. INTRODUCTION 


A. BACKGROUND AND MOTIVATION 

In the 1980's, there were tremendous advances in the 
theories that describe nonlinear systems in many fields of 
science and engineering [Ref. 1]. These nonlinear systems 
vary in complexity from a simple driven pendulum, to global 
weather systems. This thesis studies the nonlinear system of 
wind wave interaction on a "deep ocean" using an existing 
wind wave tank, Upper Ocean Simulations Facility at the 
Naval Postgraduate School (NPS). 

A theory describing the power spectrum of ocean surface 
wave crests, independent of the fetch (the distance over 
which wind blows over water) and duration of the generating 
wind, was developed by Phillips [Ref. 2]. The Phillips 
Spectrum, based on dimensional analysis, shows a 173 
spectrum for large angular frequencies, 0. 

Zakharov and Filonenko later described the stochastic 
oscillations of the surface of a liquid with a 07% power 
law spectrum [Ref. 3]. In the late 1980's, Larraza applied a 
07% power law spectrum to wind waves, flicker noise, solar 
wind, and classical second sound [Ref. 4]. 

Rene SA thirty year history of reported field and 
tank wind wave and driven gravity wave spectrum 


measurements, summarized nicely by Donelan, et. al [Ref. 5]. 


When Dr. Larraza came to the Naval Postgraduate School, we 
had an opportunity for experimentalist and thecorise 
collaborate in the laboratory while exploring his recent 
power law spectrum and second sound theories in a controlled 
environment. 

There is also a history of laboratory controlled, an 
wind wave tanks, wind driven phase speed measurements. Some 
recent experiments which have been done to understand this 
phenomenon has been either flawed [Ref. 6] by aliasing due 
to excessive spacing of probes or were not designed to look 
for the effects of a turbulent restoring force [Ref. 7, 8]. 

Wave turbulence are those processes in open driven 
systems which are dominated by a random redistribution of 
energy (generally due CO time reversible nonlinear 
processes) with effectively infinite degrees of freedom 
[Ref. 9]. Significant theoretical progress in the study of 
random wave systems driven far off equilibrium has lead 
Larraza, Garrett, and Putterman to the startling predicumen 
of a second surface wave mode [Ref. 10]. This "second sound" 
mode corresponds to propagating compressions and refractions 
of the surface roughness. 

B. SCOPE 

This is the first thesis involved in the study OE 

nonlinear wave turbulence using the existing Upper Ocean 


Simulations Facility at the Naval Postgraduate School. The 


primary purpose of this thesis is to develop, test, and 
calibrate a wave height measurement system that may be used 
in subsequent research efforts studying various wave 
turbulent effects. 

There were two wave height measurement system 
investigated in this thesis. A two-wire capacitance system 
and a four-wire conductance system were independently 
developed and tested. 

Preliminary measurements were made to verify the 
dispersion relation for driven deep water gravity waves 
[Ref. 11] 

o 

BES al 
where g 1s the acceleration of gravity and k is the wave 
number, and an investigation of the wind wave power spectral 
density and quality factor Q of the tank were conducted to 
help identify any tank limitations. As 1t stands, the second 
sound theory requires a power spectrum that decays as o>. 
All these measurements are designed to lead to successive 
experiments which may confidently measure wind wave second 
sound. 

This thesis has four chapters which describe the various 
aspects of the experiment. Chapter II discusses the theories 
of the physical properties which will be measured, and the 


theories of the various measurement techniques used. Chapter 


III describes the wind wave tank facility and it's 
capabilities, the experimental apparatus construction, 
development, and the various calibration techniques for both 
systems, and the various experimental measurement techniques 
used. Chapter IV will present measurement results. Finally, 
Chapter V will discuss the conclusions of this work, and 


suggest directions for future research efforts. 


11. THEORY 


A. INTRODUCTION 


This chapter discusses theories related to gravity 
waves, in three broad categories. These categories are 
variables, measurements, and techniques. The chapter will 
define surface gravity waves and energy density variables 
and then derive surface gravity wave dispersion 
relationships. The theories for cylindrical capacitance and 
four wire measurements will then be discussed. Finally, 
least squares fit and deviation data reduction techniques, 
Hewlett-Packard (HP) Cross correlation measurement 
techniques, phase speed measurement, wavelength, and quality 
factor measurement techniques will be discussed. 

B. SURFACE GRAVITY WAVE AND ENERGY DENSITY VARIABLES AND 


DISPERSION RELATIONS 


The [free surface of a liquid in equilibrium in a 
Gravitational field is a plane. If, under the action of 
some external perturbation, the surface 1s moved from 
its equilibrium position at some point, motion will 
occur in the liquid. This motion will be propagated over 
the whole surface in the form of waves, which are called 
gravity waves, since they are due to the action of the 
gravitational field. Gravity waves appear mainly on the 
surface of the liquid; they affect the interior also, 


but less and less at greater and greater depths. [Ref. 
TA] 


A reciprocal relation between a deep water surface 
gravity wave phase speed and frequency will be derived in 


this section. Deen water is defined as a semi-infinite fluid 


below a free surface. The fluid depth H is much greater than 


A E 
the reduced wavelength = An intuitive feel for the surface 


gravity wave dispersion relation was developed by Sir M. J. 
lever hi) 15 using simple energy arguments (Rayleigh's 
method) [Ref. 13]. A more rigorous derivation using invisid 
equations fOr incompressible flow in a constant 
gravitational field can be found in “Fluid Mechanics" by L. 
D. Landau and E. M. Lifshitz [Ref. 14]. Both developments 
will be presented here. 
1. Rayleigh's Method 

The gravitational potential energy PE per unit 
horizontal area, at any point on the surface, is found by 
integrating the potential energy per unit volume from the 
surface to the bottom. The maximum excess potential energy 
per unit horizontal area over the displacement of the 


undisturbed free surface (C€=O)is 


PE = > pat’, 
EN ZN 
where p 1s the volume density of the liquid. 
The kinetic energy KE per unit horizontal area of a 


Sinusoidal wave is 


i E 
=“ ik 
KE > pk (se) ; 


o 
where t is the time 
The mechanical energy E per unit area of the system 
is E = KE + PE. The deep water surface gravity wave 
dispersion relation is found by setting the maximum value 
of equation (2.1) equal to the maximum value of equation 
2), 
o? = gk, 
E 1 
which can be further reduced using phase speed c, c = w/k, 
and the relationship between angular frequency and natural 


frequency, ® = 2nf, to 


ace 
Df 


Gi 

EGI ia 

The above simple result assumed that there was an 

exponential decay of velocity away from the surface which 
restricts equation (1.1) to deep water. 

2. Invisid Equations 

The derivation using invisid fluid equations begins 

with Newton's Second Law, F = ma. Bold letters represent 

vectors. The equation of motion of a fluid is determined if 

this force is equated to the gradient of the pressure p on a 


ME” volume, F = ma = Vp. The equation of motion of a fluid 


was developed by L. Euler in 1778 [Ref. 15]. EU 


equation in a gravitational field is given by 


Sy 


Ót 


Vp 


+ (v: V)jv = - +g. 
p 


En Ta 
where v is the particle velocity. 

Fluid particles in a gravity wave travel on the 
order of the amplitude A of the wave during a time interval 
of period T. The particle velocity ıs of the order A/T 
velocity varies noticeably over time intervals of the order 
T and distances of the order of the wavelength A in the 


direction of propagation. The time derivative of the 


Ov. . 
velocity Se is of the order v/T, and the space derivative Vv 


is of the order v/A. The left side of equation [24005 


approximately equivalent therefore to 
Ov 


eee a 


dt 


Hp 


il 
7 + 
Eqn. 2a 
We can neglect the (v:V)v term if the amplitude of the 
oscillations in the gravity wave are small compared with the 
wavelength, A << A, and we have potential flow. 
The equation of continuity in a potential flow of an 


incompressional fluid is 


Vê4 = 0, 
Eqn. BM 


— ia 


i.e. Laplace's equation for scalar potential b, v = Vd. 
Integrating equation (2.4) after neglecting the (v-V)v term 
and a little algebraic manipulation, gives 


ôb 


p = -Pgz - P St * 


EGO | 
Let C = Č(x,y.t) be the z coordinate of a point on 
KN Na ae WIT a constant pressure pg acts on the surface, 


equation (2.7) reduces to 


0 = gl - ar 
Egn. 2.8 
The z coordinate of a point on the surface is now 
c=- -4% 
g dt z= 
Egn- 2-9 


The vertical component of the velocity of a point on 


eme surface is 


_% 
YZ = Be’ 
Ema? 10 
but vz = >. so thart 
ES 
ÓZ Zo Es 
Bem 2 IM 


Substituting equation 2.9 into ZM MEM 


5d 1820) _, 
Sz `g z= ` 


Egn. 2 
The system of equations which determine the motion 


in a gravitational field are: 


Vê) = 0, 
Eqn. T 
56 1 66) 
Oz IJ 512 > 
Eqn "Za 


Assume that the velocity decays exponentially away 
from the surface at a rate k. The velocity potential of the 
deep water surface gravity wave propagating in the x- 


direction and uniform in the y-direction is, 


db = A et Zcos (kx-0t). 


Eqn. 24 
Substituting equation 2.13 into the boundary condition oaei 


equation 2.12, we obtain 


Eqn. Zee 


TO 





or 
02 = GKF 

ESmes 1.1 

This gives the same deep water surface gravity wave 

dispersion relation between wave number and angular 
frequency, as developed using Rayleigh's method. 

The deep water approximation will start breaking 

down when the bottom of the tank interferes with the 

Sedution given by equation (2.13). This will occur if the 


¿KZ 


term 1S non-negligible at z -H, H being the depth of 


1 
the water, which occurs when k < H Eor our tank, H = Im so 


the deep water approximation 1s good for 


on; al 


A >> + = 21H = 6.3 m and Ji NS EO Hz: 


C. TRANSDUCER PRINCIPLES 

Understanding the electrical properties that are 
proportional to changes in surface water height are an 
essential part of this thesis. The electrical properties and 
measurement techniques used in this work will be discussed 
next. 

fewer electrical circuit, complex electrical impedance 
1S proportional to voltage changes when a constant current 
ES applied to the circuit. The first type of transducer 
described in this thesis is based on applying a constant 


voltage and measuring the change of capacitance with respect 


L1 


to the length of a cylindrical capacitor as a funct tone 
water height. 

The second type of transducer described in this thesis 
1s based on applying a constant current and measuring the 
change in the magnitude of the complex impedance with 
respect to an effective cross sectional area of a volume of 
water under test. The impedance is inversely proportional to 
the water height. 

This section will detail both transduction properties 
and the associated measurement techniques. 

1. Cylindrical Capacitance 

A vertical insulated wire penetrates the surface of 
water approximately forming a cylindrical coaxial capacitor. 
The wire is the inner conductor, the insulation is RE 
dielectric, and the water bath acts as an outer conductor of 
changing length. 

To ensure that our system is behaving as expected, 
the theoretical limiting capacitance of an ideal cylindrical 
conductor must be compared to the measured capacitance of 
the cylindrical probe wires. The capacitance of an ideal 


cylindrical conductor is given by ARSS 


2 
c = 1, 
1n(2) 
a 


Ban. 2 


where 


12 


Se permitrivitey,.e = Keo; 

K is the relative dielectric constant of a material; 

E> 1S the permittivity of free space; 

Wwe cme length of the conductor, proportional to the 

height of the water; 

In is the natural logarithmic function; 

bis the radius of the outer charged surface; and 

ais the radius of the inner charged surface. 
Mme permittivity of free space is £€ = 8.854 x on! 
farad/meter. 

The resistivity of the water in our tank is 
Surprisingly low. A few inches of 0.2 mm diameter wire has a 
resistance to infinity of only 20 ohms. The capacitive 
impedance of our probe is many times higher, so that the 
water can be taken to be a perfect conductor. Electrical 
contact to the water is made in a two-wire measurement 
configuration through a second submerged wire at "infinity", 
many probe wire diameters away. 


A measurement of the capacitance gives the height, 
. mz, ÔC | ee 
with the theoretical sensitivity <7 for the final capacitive 


ôl 


| 2TE€ 
test probes given by m~ A E / 10 
1n(2) 


a 


The two-wire capacitive scheme suffers from an 
unknown, unsteady, artifactual complex impedance between the 


probe and the bath at the very small length scales of the 


iV) 


order of a probe diameter. For example, the water meniscus 
slips and sticks to the probe because of a surface film that 
collects on the probes, instead of sliding freely (SAS 
water height oscillates. This film can itself add an unknown 
series resistance. And it is also possible that the 
conductivity of the water in the immediate vicinity of the 
probes is changing because DÉ poorly understood 
electrochemistry at the metal-water-ion interface. For these 
reasons, we developed a four wire probe, based on the four 
wire impedance measuring technique commonly used in 
electronics and solid state physics. 
2. Measurement Technique Models 
The physical electrical connections between an 
analyzer and a device under test (DUT), water in our case, 
using a two-wire measurement technique can be modeled as 
unknown variable resistors in an equivalent electrical 
circuit as shown in Figure (2.1). A constant current 
measured and applied to the DUT and the voltage v is 
measured across the series combination of the unknown lead 


impedance and the DUT impedance. 
<> 


Figure 2.1. Two wire measurement equivalent circuit. 
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lour wire Keng or cermmnalpair. configuration, 
measurements provide significant measuring advantage over a 
two wire measurement, for some situations. The HP-4192A and 
HP-4194A Impedance Analyzers are designed to make four wire 
measurements. The combination of equipment and measurement 
technique provides an accurate and easy measurement. 

A four wire measurement uses four connectors: High 


Current (Heur). High Potential (Hpot). Low Potential (Lpot). 


and Low Current (Leur) [Ref. 17]. The unknown variable 
impedances of the potential leads, as shown in the 
equivalent electrical GEC TE in Figure Wa are 


unimportant because there is no current in that portion of 
the circuit due to the high, nearly infinite, impedance of 
the voltmeter. Likewise, the unknown impedance of the 
Current leads are unimportant because the current flowing 
through them is known, and all this current flows through 


the DUT. 


H $7 de 
Cur Cur 


Figure 2.2. Four wire measurement equivalent circuit. 
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3. Four-Wire Conductance Technique 
Conductance 1s inversely proportional to TENS 
magnitude On the complex impedance. Conductance iS 
proportional to changes in the water height. The conductance 
is obtained using the four-wire measurement technique by 
inverting the voltage difference in the water when measuring 
the complex impedance. 
The complex impedance Z of a device under test DUT 
is given by 
“hn Were i 
Ean . "BS 
where 
p is the complex resistivity of the conductor (water in 
our case); 
l is the length of the conductor, proportional” Comma. 
separation between the potential leads in the water; 
h is the height of the conductor on the Ire 
Weff is the effective width of the conduction path 
between the wires; 
AV = Hpot a Loot i 
i is the current applied to there olt 
Figure (2.3) illustrates a four-wire resistance measurement 


in water. 
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Equipotential H H 
Surfaces 


—Wire Conductor 


pa TOS EE 
Water air - 


interface 





Figure 2.3. Four-wire resistance measurement. 


Because it is based on the four-wire technique, our 


probe is insensitive to possible conductance changes of the 


water or to other contact impedances between the wire and 


the water on the "micro scale" illustrated in Figure (2.4). 


However, the four wire technique would be sensitive to 


changes in conductance on the "macro scale" between the 


probe wires. We do not know if this is a problem, or not. 


Macroscale 








Probes 


ale 
OO O 


Microscale 


Figure 2.4. Macro/micro scales about a four-wire sensor. 
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D. LEAST SQUARES FIT AND DEVIATION ANALYSIS TECHNIQUES 

The method of least squares for fitting data to a 
straight line can be accomplished on "most handamam 
scientific calculators today. The calculator automatas 
accumulates the necessary sums to provide slope, intercept, 
and linear correlation coefficient. Unfortunately, the 
correlation coefficient is not directly related COMAS 
uncertainty in the least squares values of the slope and 
intercept, which are the quantities of interest. [For a full 
development of the method of least squares fit, see 
reference 18). 

The correlation coefficient r indicates whether there is 
causal relationship between the variables x and y, and 
therefore whether a linear fit to the data is justified. The 


correlation coefficient is defined by Taylor [Ref. 19] in 


terms of the covariance, o and standard deviations O, and 


xy 


O,- A more intuitive picture of the correlation can be 


developed by determining the slope m of a line of the data 
set when it is plotted y vs. x, and then when it is plotted 
x vs. y (which we will call m'). If there were no errors in 
the measurements and the measurement x and y obey a linear 


relation, then these two slopes would be a negative inverse 


of one another (i.e. m = -l/m') and r2 = 1, where TM 
m 

defined as r = EE If there is noise in the data, them r2 < 

ie 


18 


fio quantitative Seagmaticance of r is found when r is 
expressed in terms of probability. For any defined ro, 
Py(Irlzr,) is the probability that N measurements of two 
uncorrelated variables would give a coefficient r as large 
as ro. The probability that N measurements were derived from 
a random process involves the integration of a two 
dimensional Gaussian distribution [Ref. 20]. Generally it is 
more convenient to use a table which contains a value for r 
which would have to be exceeded for there to be a 
probability less than Pytlrbro) that the N measurements were 
generated by a random sample of observations taken from 
uncorrelated parent population. Two such tables are included 
Tay lor [Ref. 21]. 

For physical measurements, it is the error in the least 
Square values of the slope and intercept which are of 


interest. Taylor gives expressions for the standard 


deviation in y, O error in the slope, o and intercept, 


yee 
Oy. in terms of the accumulated sums which are available in 


m’ 


memory of most scientific calculators. I have reproduced 


them below: 


Oy 





Bem 2.16 


o 


2 
Sm y 2 2 Ega => 22 (xXx Yp)12xp2Yp + NI Laa 
n 





= Ly = 
n 2 2 ; 
Ei N2x,- (2x,) 

Eqn. Mm 

Op = 
Fan. WENG 

Om = 
Eqn. 200 


where n 1S a summation dummy index. 

Data sets were analyzed and errors in the fit parameters 
were determined using this technique throughout the thesis. 
Program code to determine the standard deviations in 
intercept, slope, and y using a HP-15C calculator is 
documented in Appendix A. The data was also analyzed using a 
Microsoft™ Excel file using equations (2.16) CARO 
(2 200: 

It 1s also important to determine 1f there are random or 
systematic errors when examining experimental measurements. 
Random errors are revealed by repeating the measurements. 
Systematic errors cannot be revealed this way. Deviations in 


the least-square fit data are computed by subtracting the 


20 


regression calculated dependent value ya, from the measured 
dependent value ya [Ref. 22], 
Penang = yo = Yc- 
BC 2 20 
Plotting the deviations helps reveal what type of errors are 
showing up in the measurement data. 
E. HEWLETT-PACKARD ANALYZER CROSS CORRELATION MEASUREMENTS 

"Correlation is a measure of similarity between two 
quantities" [Ref. 23]. This measurement technique allows you 
to measure a very small signal in the presence of noise. The 
dynamic signal analyzer can make this type of measurement 
well for very low frequency or single shot measurements, and 
provides a signal output in the time domain. 

There are two types of correlations: the auto- 
correlation function and the cross-correlation function. The 
results of the different correlations can be the same; 
however, the different correlations are applied to different 
types of systems. The auto-correlation function is a single 
channel measurement and the cross-correlation function is a 
two-channel measurement. 

The cross-correlation function compares two different 
Signals as a function of the time shift between them [Ref. 


24), and is given by 
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Son 


ii 
xy ET x (tC) y TCCs 


O Cam 


Eqn. 7a 
where x is one input, y ıs another input, © is Imo mou 
and t is the time shift between the inputs. 

When the same signal is present in both channels in the 
presence of noise, the noise can be greatly removed using 
this technique. The correlation analysis is made when the 
dynamic signal analyzer is making a power spectral density 
measurement: 

Power spectral density measurements are made by the 
dynamıc signal analyzer for each channel. A cross power 
spectral density is then computed internally, which is the 
product of the Fourier Transform of the two signals =i me 


25], and 1s given by 


Eqn. 29a 
where f is the natural frequency and * indicates the complex 
conjugate of the function. 

The HP cross-correlation is then computed indivegeuas 


using the cross power spectrum [Ref. 24], and is defined as 


LES 
Ra [Ga AA 


Eqn. 22a 
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F. GRAVITY WAVE PHASE SPEED AND WAVELENGTH MEASUREMENTS 

The phase speed of a sinusoidally driven gravity wave 
can be measured in a two-point measurement system using the 
output of the signal analyzer cross-correlation function and 
the distance between the sensor probes. This technique was 
used as a check of the capacitance probes. 


The phase speed of a wave is given by 


A, 
Eqn 2.24 
and the theory of deep water gravity waves gives this speed 


as 


a 


C = on’ 


Pons, 245 
Experimentally, the phase speed of a wave in the x-direction 
can be determined from the period of time T that a 
particular phase of the wave, such as a crest or trough, 
takes to move a given distance d. The phase speed is given 


by 


O 
H 
Ha 


Ban. 2.25 


where the distance d is the distance between the probes. 


Zo 


The cross-correlation output of a sinusoidal signal will 
also be a sinusoid, in the time domain. The period of a full 
or half cycle of a driven surface gravity wave can be 
determined from the cross-correlation measurement output. 
The period is determined using the method of least squares 
fit to the measured cross-correlation minimum or maximum 
data points. 

The phase speed is then a simple computation using 
equation (2.25). The driven surface gravity wave wavelength 
1s determined using the measured period and computed phase 
speed. The natural frequency is the inverse of period, i.e. 
j= L Substitute this into equation (2.24) and solve 


T 


wavelength to get 


Eqn. 2m 

G. FREE DECAY QUALITY FACTOR, Q MEASUREMENT 
The amplitude of a damped oscillator decreases 
exponentially with time. The amplitude may be represented by 


Aes) = mee 


Eqn. 2a 
where Ag is the initial amplitude and ß is called 


temporal absorption coefficient. The length of time oa 
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damped oscillator to decay to 1/e of the initial amplitude, 


is called the decay time t, defined by 


TH lh 
| 
E 


cem 2.28 
where m is the mass and Rm is the mechanical resistance 
Peer. 26). 

Substitute equation (2.28) into equation (2.27) and take 


the natural logarithm of the equation to get 


i 
IAS fe) ie E 
Eqn. 2.29 
Using a change of variables, letting y = ln [A(t)], m = =. 
and b = ln [Aj] we get the equation of a line, 
y =mt + b. 
Bam. 230 


For a continuous frequency system, such as the spectrum 
of a wind wave tank, the amplitude vs. time can be measured. 
Plotting ln [A(t)] verses time, b and m can be determined 
using the method of least squares fit. The decay time can 


then be computed for any frequency fo by 


EGQite 2 eo 
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Another method of reporting this decay is in terms of 


the quality factor, or "O of the system". The quality A 


of a damped oscillator is defined [{Ref. 27], 


Ean 


where œ> and 607 are average half power frequencies, above 


and below resonance Oo, respectively. 
I will not discuss the physics of a damped oscillator. 


This is discussed in "Fundamentals of Acoustics" [Ref aan 


The Q can be expressed in terms of the mechanical constants 


of the system, given by 


Ean 2m 


Substituting the decay time from equation (2.28) gives 


Eqn. Ji 


Note that Q is a dimensionless quantity. 
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III. EXPERIMENT 


A. INTRODUCTION 

The goal of our experimental program, of which this 
thesis is the first part, as to measure second-order 
(nonlinear) effects on deep water surface gravity waves. A 
very sensitive and linear measurement system 1S required so 
that transductTon antmhacts, such as non  linearity, 
hysteresis, noise, etc., will be minimized. Interpretation 
of second-order effects will be much more reliable 1f we can 
ensure that we are using a linear, low noise, high 
sensitivity, wave height transducer system. 

There were two wave height transducer systems 
investigated in this thesis. The two-wire capacitance system 
was developed and tested until it was determined that the 
system was unsuitable for these applications. A four-wire 
conductance system was then developed and tested. The 
transducer systems development and testing are presented in 
parallel. 

The experimental apparatus can be divided into four 
major components, as shown in Figure (3.1). The first 
component 1S the wind-wave generating tank system. The tank 
system was originally built in the 1960's by Professor 
Herman Medwin for use in experiments to measure sound 


scattering produced by wind driven waves on rough surfaces 


Za 


[Ref. 29]. It was upgraded for these experiments. The second 
component is the sensor. Two sensor types were tested in 
these experiments. The third component is the associa 
electronics designed and built for each sensor. The final 
component 1s the off the shelf, general laboratory 


instrumentat iong 


Wind-Wave Analyzing 
Sensor . 
Tank System Equipment 


Figure 3.1. Experimental block diagram. 





Associated 
Electronics 





This chapter will describe the general calibration 
techniques used, some of the design criteria and detailed 
construction of the experimental apparatus components (1.e., 
sensor, associated electronics, and the wind-wave tank 
system), and the specific calibration of all the apparatus: 
Finally, the techniques and experimental procedures used for 
the deep water surface gravity wave measurements will be 
detailed. 

The thesis will have several charts and tables that are 
referenced to the thesis lab books. The system will be lab 
book number, page number, represented by (B#P#). 

B. WIND-WAVE TANK SYSTEM 
1. Tank Facility 
This section will” descaie- oe. capabilities and 


dimensions of the existing wind-wave tank and upgrades, the 


228 


first experimental component. The wind-wave tunnel section 
MS feet long, búilt primarily above ground, as shown in 
megure (3.2). The tunnel is constructed of 3/4 inch plywood, 
meer glassed for waterproofing, swith 19 (2' x 4' x 0.5")(1 
x w x t) plate glass view ports. The water cross section of 
the tunnel is (4 x 3.5 feet) (w x d). The wind cross section 
em the tunnel is (4 x 0.5 feet) (w x d). The depths are, of 
course, only approximate, since they depend on the quantity 


of water placed in the tank. 
Top View 


| a ez! a kdi a 
as! EE 


10 ' 
a | Beale 


Vievports 







ee Housing Anechoic 
APRA A PA t 


Front Side View 


Figure 3.2. Wind-wave tank. 

The anechoic section is (10 x 10 x 9.5 feet)(1l x w x 
NITIK into the ground, also shown in Figure (3.2). This 
Æ tion of the tank is constructed of six inch reinforced 
concrete and has (4" x 4" x 9.5')(l xw xt) redwood wedges, 
used on the four sides and the bottom to absorb underwater 
sound in previous experiments. The top of this section has 
two removable 3/4 inch reinforced plywood covers for 


personnel safety. 
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2. Small Wave Guide 
A paddle driven at large amplitude, can generate 
parametric waves that run transverse to the directly driven 
wave. These waves interfered with the dynamic calibration of 
the four wire probe. You can cutoff these unwanted waves 
using a wave guide. A small wave guide was constructed and 


placed in the main tank as shown in Figure (3.3). 





Plexiglas— 


Figure 3.3. Small wave guide. 

This wave guide has only two sides and a small 
paddle at one end. The top, bottom, and end opposite the 
paddle are open to the main tank. The minimum width of the 
wave guide was determined by the minimum width of the 
conductive probe, which was calibrated in this wave guide. 
The cutoff for this wave guide Ts Bos a 

One side of this wave guide is a (2" x 6" x 10')(t x 
w x l) pine board. The other side is a (0.060" x 2' x 
x w x 1) Plexiglas plate. The Plexiglas is used so that the 


wave structure in the wave guide can be seen. The two sides 


3 
are held together with six, eight inch long, g16 bolts, 


i ; 3 
five, four inch Teflon spacers between the sides, and s1x, 8 


-16 hex nuts. 
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AA AEREA x |) pine board” was used 


Bor the paddle. An eight inch long, 2-16 bolt went through a 


half inch hole that was drilled at one end of the paddle and 
two (5 x 5.5 x 16 inch)(t xw x 1) pine boards that were 
also bolted to the wave guide sides. A wave driver is 
connected to a 16 inch L-channel which was bolted to the top 


of the paddle. 


3. Wind Generation 


Aju 


The wind —source  consist of five pairs of 


horsepower centrifugal fans on the top of a (39 x 48 x 41 
inch) (1 x w x h) housing. This housing 1s attached to the 
tunnel between locations one and two in Figure (3.2). There 
are three power switches for the fan bank resulting in seven 
combinations of powered fans. 

The wind was measured approximately four feet from 
the wind-water interface, approximately in the center of the 
wind tunnel cross sectional area, at location two in Figure 
(3.2). Wind speed was measured using a hand held Turbo- 
meter™ manufactured by Davis Instruments of San Leandro, 
CA. The uncalibrated, measured wind speeds are reported in 


Table (3.1). 


orl 


TABLE 3.1. WIND AND SURFACE DRAG SPEED ATALOCATTONHHEN 


(BAP IA). 


Surface 
Fan Drag Speed 
Combinations (m/s + 0.1) (m/s + 0.01) 
Fi Too small to 


measure 
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4. Wave Driver 
The wave driver 1s a gravity wave source. Three wave 
drivers will be described in this section. 
a. Medwin Wave Driver 
The Medwin wave driver was developed and used in 
previous experiments in the tank system. The Medwin wave 
driver has a fixed amplitude and four drive frequencies. It 
consists of a forcing wedge and belt driven electro- 
mechanical plunging system. New wave drivers were developed 
for these experiments after initial tests 1n the tank 


determined that the Medwin wave driver was nonlinear itself. 


ay 


The forcing wedge was a right triangle wedge 
that is pushed and pulled up and down in the water, Figure 
WA). “INEA Io ue Inch right triangle wedge is 
constructed using = inch plywood sides, braced with à inch 
plywood on both ends and in the center, and open on the top. 
There are one and a half inch guide strips on the outside of 
both ends of the wedge that are inserted into tracks built 
in the end of the tank just behind the fan housing. The 
tracks help direct the wedge movement up and down into the 
water. The entire wedge 1s treated with waterproof paint for 
submersion into water. Sand is uniformly placed inside the 
wedge to increase the weight in an attempt to increase the 
wedge stability. 

The plunging system consist of an electric 
motor, a fixed mechanical speed reducer, fly wheel, and two 
plunging rod assemblies shown in Figure (3.5). The electric 
motor 1S a Dayton Electric MFG. CO., model number 5K452C, 
half horsepower, 1725 RPM, 60 Hz, single phase, capacitive 
start motor. The motor is mounted on standard (1 x 1.5 x 10 


mem) motor mounts. The motor mounts are attached to a (11 x 


OO [to 


8 x 3 inch)(1 x wx t) aluminum plate using four, one inch 


Reno, =-16 Ps plate is attached to a (13 x 22 x = 


Ren) (1 x w x t) aluminum bottom plate using an eight inch 
piano hinge and a four and a half inch long, 5-13 ¡NAAA dis 
construction allowed the position of the motor to be easily 


changed for belt adjustment and/or replacement. The large 
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bottom plate is shock mount to two (3.5 x 3.5 x 2475 amenas 
x w x 1) redwood boards. These boards are clamped to 
aluminum L-beams that span the top of the tank, just behind 
the fan housing at location one in Figure (3,2) The 


plunging system is mounted above the water directly above 


Y, 7 “ 0 Y 
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Figure 3.4. Medwin wave plunger. 
The mechanical speed reducer is an Alling-Lander 


Co., Inc., model number CU-20, 20 to 1 ratio reducer An 


aluminum ten inch diameter flywheel is attached to the 2 


i ; 144 
inch diameter output drive axle. An eight inch long, g inch 


diameter stainless steel rod, with two and a half inches of 


5-13 threads is connected to the output flywheel. 
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Figure 3.5. Medwin drive system. 


A Craftsman™, catalog number 9-1631, ( So 


ire 


inch) engine belt is used to connect the electric motor to 
the mechanical speed reducer. The speed reducer is mounted 
to the bottom plate also. The belt is used on two, four 
position, flywheels. Positioning the belt in any of four 
flywheel combinations determined which of the four output 
frequencies the system generated. 

The plunger rods assembly consisted of two (2 x 


2.5 x 1 inch) aluminum bearing housing, which are attached 


ss: 4 
to the 3 inch diameter rod on the flywheel. The housings 


have a 2-16 threaded hole that a 14.5 inch threaded rod may 


be screwed into and held with a lock nut at the housing. A 


seven inch turn buckle is attached to the other end of the 


> 


threaded rod. A two and a half inch eye bolt is attached to 
the other end of the turn buckle. A Six Incas 20 bolt 
with nuts goes through the eye bolt and the center pieces of 
plywood in the wedge, completing the plunger system. 

The non linearity was not quantified, but was 
readily apparent when watching the erratic plunging action 
and resulting waves of this driver. The wedge would not move 
up and down in a uniform motion. It erratically Got Cama 
in the guides due to the fact that the force from the 
plunger rods was not always normal to the wedge itself. This 
caused the wedge to be pushed or pulled from side to side as 
well as up and down. 

b. Yarber Wave Driver 

The Yarber wave driver iS continuously variable 

in both frequency and amplitude. It has a forcing paddle and 


a direct drive electro-mechanical plunging system. 


The paddle is a (12 x 44 x = inch) (w x | = 


pine board with a (15x 4 2 inch) tab, which is Just ree 


of the center line as shown in Figure (3.6). The paddle is 
attached to two (2" x 4" x 4')(t x w x 1) boards witha 
three and a half inch hinges. The two by fours are clamped 
to the end of the wave tank, using two six inch C-clamps. 
The paddle system consist of an electric motor, 
a direct drive speed reducer, and a single plunging rod 


assembly shown in Figure (3.7). The motor is a General 
Electric Co., model number 5KC45PG1140H, 2 horsepower, 1725 
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RPM, 60 Hz, single phase, capacitive start motor. The motor 


is connected to a Zero-Max® direct drive, model number OXl, 
CCW output, zero to 400 RPM continuously variable speed 
mechanical reducer. Full specifications are found in 


Reference 30. 





Figure 3.6. Yarber wave plunger. 


The ten inch flywheel and plunger’ rod 
assembly were removed from Medwin's wave driver and 
connected to this system. The plunger rods are attached to 
the flywheel using a new (12 x 2 = A ASE ct) 
slotted aluminum adapter. The slotted adapter provides a 
continuously variable adjustment of the driving amplitude. 


The plunging system was mounted to a one cubic foot wooden 


E 


box, that was clamped to a cement laboratory table. The 


cement table was placed as close to the tank as possible. 
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Figure 3.7. Yarber@@rive syste 

The Yarber wave driver would change frequency at 
a rate of approximately 0.002 Hz per minute (B2P61) when the 
system ran continuously for over two minutes due to 
overheating in the mechanical speed adjuster. The amplitude 
and frequency were not adjustable during a measurement. 
This combination of things caused us to develop ¿MMS 
drive system. 

C. Keolian Wave Driver 

The Keolian wave driver is also continuously 

variable in both frequency and amplitude. It uses the same 
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forcing paddle assembly as the Yarber wave driver and an 
electric plunging system. 

The plunging system consist of a function 
generator, power amplifier, electric shaker, and a single 
plunging rod assembly shown in Figure (3.8). A HP-3314A 
function generator is connected to an Acoustic Power 
Systems, Inc. APS-114, dual-mode power amplifier. The 
amplifier is operated in the current mode. The amplifier is 
connected to an APS-113 Electro-SELS® shaker. A Kikusul, 
COS-5060A, 100 MHz oscilloscope 1s used to monitor the 


amplifier input and output. 


HP-3314A COS-5060A 





© + è ù q 






Figure 3.8. Keolian drive system. 


The plunging rod assembly consist of a 14.5 


i 3 A 8 
Weh, g 16 threaded rod with a one inch ball joint on each 


end. One ball joint was connected to a three inch aluminum 


L-bracket, which was bolted to the shaker using 5 GD Long; 


Í de 
4720 bolts. The other ball joint was connected to a L- 


bracket on the paddle. 
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The wave amplitude and frequency are easily 
adjusted at any time using the synthesizer/function 
generator. Complex wave forms and noise can also be applied 
to the paddle using the amplifier and shaker. This will be 
important in latter experiments. 

5. Surface Drag Measurements 

Surface drag measurements were made at various wind 
speeds, approximately four feet from the wind insertion 
point, at location two in Figure nam 

A special apparatus was constructed to hold the 
wires and associated electronics of the capacitive probe 
setup in the tank. We call this apparatus the Harp. The Harp 
is constructed of two 47.5 inch long aluminum U-channels, 
connected by two 30 inch long, 1/2 inch diameter aluminum 
rods. The upper U-channel is 2.25 inches wide by one inch 
high. The lower U-channel is 0.9 inches wide by 0.6 inches 
high. The wires were mounted using number one rubber bands 
attached to one end of the probe and the lower piece of the 
Harp. The upper ends were attached to a guitar adjuster, 
which was mounted to a 1/4 inch Plexiglas strip that was 
connected to the Harp. An electronics box was screwed onto 
the Harp. 

The Harp was put into the water with a single 
Capacitance probe and a clear plastic ruler attached to the 
Harp with a three inch C-clamp. The ruler and probe were 
separated by ten centimeters at the air-water interface and 
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the line formed was normal to the wind direction. Lint was 
placed on the water upwind of the ruler. The dust travel 
time between the two points was measured using a Casio Alarm 
Chrono G-Shock digital watch. The measured surface drag 
speeds are reported in Table (3.1). 
C. CALIBRATION TECHNIQUES 

1. Introduction 

There are two types of calibration techniques which 
were performed in this thesis. The apparatus and procedures 
for two calibrations will be described in this section. 

The first type of calibration was a static 
calibration which was used primarily during the probe 
selection. Static calibrations were also conducted to 
establish system baseline sensitivities. 

The second type of calibration was a dynamic 
calibration. The dynamic calibrations were the most 
important for the thesis because they set the sensitivity 
limits for the sensor system for all further measurements. 

2. Static Calibration 

The static calibration technique had three component 
parts as shown in Figure (3.9). The specific parts used as 
the water or electronics components could be changed, 
depending on the specific test conducted. A Starrett Co., 
No. 254, 26 inch master vernier height gage was used in all 


ceat nic calibrations. 
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The static calibration technique is simple. The 
sensor or sensor system was somehow attached to the height 
gage. The necessary measurement electronics were connected 
to the sensor or sensor system. The sensors were immersed in 
water. The sensor parameters or sensor system electrical 
outputs were measured, verses the measured relative 
immersion determined uSing the height gage. A sensitivity 


was then computed, in the form 


A Our pu 
A height 


Eqn SEM 






Figure 3.9. Static calibration components. 
3. Dynamic Calibration 

The dynamic calibration techniques also have three 
component parts as shown in Figure (3.10). All of these 
components could be changed depending on the specific test 
conducted, and components could be combined during certain 
tests. The dynamic calibrations were limited to sensor 
system tests. 

There were two movement possibilities for the 


dynamic calibrations. The sensor system could be moved in 
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the water or the water could be moved past the sensor 


system. Dynamic calibrations for both movements were 


attempted. 






Figure 3.10. Dynamic calibration components. 

A subjective decision was made that the best dynamic 
calibration was an dynamic calibration for a constant 
amplitude traveling water wave, moving past the sensor 
system. This decision was made because the traveling wave is 
a gravity wave. Gravity waves and wind-water waves were what 
the sensor system was designed to measure. The other dynamic 
calibration techniques attempted are described fully in 
Appendix B. 

The sensor system was placed in the small wave guide 
describe in chapter III-B-2-b. A static calibration was 
conducted to establish the baseline sensitivity. A 12 inch 
Mena, 0.03] inch diameter brass rod was scotch taped to a 
Barrett Co., No. A LH, three axis micrometer. The 
micrometer was bolted to the main tank approximately three 
feet from the paddle and directly above the small wave 


guide. 
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The brass rod was raised 0.354 cm above the 
unperturbed water. This would provide a visual point with a 
0.25 cm RMS maximum for the traveling waves. 

The Keolian wave driver was attached to the paddle 
of the small wave guide. The drive amplitude was adjusted so 
that the traveling wave amplitude was constant, + 0.1 cm, as 
determined by the peaks of the wave just touching the bottom 
of the brass rod for several frequencies below 4.0 Hz. The 
small wave guide cutoff for the paddle generated parametric 
waves was 3.95 Hz which limited the frequency of the 
calibration. The sensor system electrical output for a 
constant amplitude traveling wave was measured, verses 
applied drive frequency. A sensitivity was then computed, in 


the form 
A output 


A frequency 
Ean e 


D. WAVE HEIGHT PROBES 

The second experimental component is the wave height 
sensor probe that we need to make the surface gravity wave 
measurements. The goal is to develop a height sensor that is 
linear, very sensitive, and noise immune. A variety of wave 
probes have been used by experimenters using changes in 
Capacitance and resistance [Ref. 5, 7]. Two sensor probes 
were tested in this investigation. The constraints in the 
probe selection, test procedures, results, and final Rana 


selection will be discussed in this section. 


ad 


1. Capacitance Probe 

The capacitance probe investigated in this thesis 
utilized the two-wire measurement technique described 
earlier. The variable capacitance measured is determined by 
the probe immersion into water. The capacitance probes were 
all labeled and marked T-(#) so that repeated measurements 
could be made for analysis for each probe. The probes were 
placed on the Harp for all capacitance test conducted in the 
wind-wave tank. 

There were two major categories of constraints on 
the capacitance probe selection. Optimizing the various 
electrical and hydrodynamic constraints was an interesting 
challenge. 

There are three electrical constraints on the probe 
when it is in the water. The first constraint is that the 
resistance between the probe and water should be nearly 
Er inite. 

The second constraint is that the capacitance per 
unit length should be as high as possible to maximize the 
sensor sensitivity. From equation (2.15), the capacitance 
increases as the radius of the outer charged surface 
approaches the radius of the inner charged surface. We 
wanted as thin a dielectric as possible. Meeting this 
constraint would establish an “ideal" electrical capacitive 


Component. 


45 


The third constraint is that the capacitance 
sensitivity, should be constant over the length of the 
cylinder, to remove unwanted systematic errors, HM 
required a uniform dielectric. 

There were also two interrelated hydrodynamic 
constraints on the probe that also require optimization. The 
first constraint was that the probe drag needed to be 
minimized to reduce the induced turbulence in the flow 
field, caused by the probes. This also required a thin 
conductor and dielectric. Second, the effect of 
wetting needed to be minimized so that the change of 
Capacitance closely followed the water height, reducing 
hysterisis. 

The optimization of the hydrodynamic constraints had 
to be accomplished on two time scales. The period of a 
forced surface gravity wave is the first time frame. This 
time frame is important, as it sets the low frequency limit 
of the sensor. The period of an experimental data run is the 
second time frame. The capacitance sensitivity, must remain 
constant during this time frame to prevent systematic errors 
in the measurements. 

Finding the right sensor seemed like it would be an 
easy job. First we needed to find an electrical conductor 
that had nearly infinite resistance, that could be used as a 


conductor/dielectric/water capacitor. 
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Determining what wire should be used was the first 
order of business. The first test was to find a linear and 
highly sensitive probe. The static calibration technique 
discussed in chapter III-C-2 was used to measure the probe 
sensitivity. The sensor probes were submerged in water taken 
from the tank and contained ina (2.25 x 17.5 inch) (diameter 
x h), 1000 ml, Pyrex™ graduated cylinder. The capacitance 
of a submerged probe was measured using a HP 4192A LF 
Impedance Analyzer. 

Measurements were made with the HP-4192A using a 
30.0 kHz, 1.0 volt drive signal. The 30.0 kHz signal is used 
because comparable frequencies will be used in latter 
experiments. A 5.0 ounce plum bob is used to straighten the 
probe and insure that it 1S perpendicular to the air-water 
interface. Figure (3.11) is a diagram of the experimental 


apparatus. 


HP-4192A 










starrett Co. 
Master Vernier 
Height Gage 
No. 254 


_Pyrex 1000 ml 
Graduated Cylinder 


Plumb Bob 
Tank Water 


Figure 3.11. Static capacitance sensitivity experimental 


setup. 
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The first capacitance probes that were tested were 
all made from magnet wire. Several coatings were tested on 
the magnet wire. There were eight probes built, using 22 
gauge magnet wire manufactured by Belden and Consolidated. 
The ends of half the probes were dipped in blue Testors 
model paint. The ends of the other half of the probes were 
coated in BIPAX, TRA-BOND, BA-2106, 5-minute epoxy. Half of 
the probes were coated in Bulls Eye Shellac, while the 
remaining probes only had their manufacturer's insulating 
coating. Neither the capacitance nor the dissipation were 
constant, at a constant immersion depth, for a length of 
time comparable to an expected gravity wave measurement. 

A second capacitance probe was built, but not fully 
tested because the measurement technique was changed. A 33 
inch long, 0 loan ca z inch OD. glass capillary 
tube formed the dielectric for this probe. A Malco S-50 
Micro dot connector, P.N. 031-0001-0001, was soldered to 22 
gauge tin copper wire (Aztec Wire & Cable). One end of the 
glass was melted closed and the wire was place in the 
Capillary tube. The connector end of the wire was attached 
to the glass with several layers of shrink tubing. 

The final wire tested and selected for use was Omega 
Engineering, Inc., TF-010, 10 mil diameter, Teflon coated 
copper thermocouple wire. The Teflon has a dielectric 
constant K = 2.0 to 2.1 [Ref. 31]. It was selected for the 
sensor probe because its capacitance was constant over time. 
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Figure (3.12) shows measurement samples of capacitance 
verses time for shellacked magnet wire (T-2) and Omega 


thermocouple wire (T-15). The shellacked magnet wire 


4.5 PE The 


UE Omega 


ES 
min’ 


capacitance changed at a rate of 


thermocouple wire capacitance changed at a rate of 0.013 


* Shedeched Megnet Wire . 
O Omega TF-101 Thermocouple Wire 


0.1 


Delta C/Co 0.06 


0.02 





Time (min) 


Figure 3.12. Capacitance of immersed probe wire vs. Time 
(EBEG). 

The immersed end of the Teflon wire was looped and 
tied to the plum bob or other tensioning device. The knot 
was covered with Polymeric Systems, Inc., PSI 601 Clear, 
100% Silicone, RTV Adhesive Sealant. Approximately > inch of 
electrical heat shrink tubing was shrunk over the knot and 
RTV. The "dry" end of the probe was tinned and soldered into 


a male, 0.062 inch, terminal (Digi-Key part number WM1000). 
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Test were performed in tank water and in mercury 
using the static equipment setup shown in Figure (3.4). The 
mercury was held in a different container for safety 
reasons. The mercury test should give results as close to 
theory as experimentally possible. 

The theoretical capacitance per unit length for this 
geometry using equation (2.15) is Ctheory = 2.40 + 0.01 
pF/cm (B2P73). The measured capacitance sensitivity was 2.2 
+ 0.1 pF/cm in mercury and 2.4 + 0.1 pF/cm 1n tank War 
shown in Figure (3.13), and was linear as a function or 


immersion in Ppecth cases. 
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Figure 3.13 Capacitance of Teflon wire in different 
media vs. Relative immersion (B2P22/102). 
A wire capacitance static calibration was performed in 


the wind-wave tanks using precision blocks for the height 
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gage. The blocks were 1.7 + 0.1 cm high. The measured 
@apacitance sensitivity for this test was 2.2 + 0.1 pF/cm, 
as shown in Figure (3.14), and was linear as a function of 
immersion. 

A final static capacitance sensitivity test was 
performed in the large wind-wave tank using a HP4194-A 
Impedance/Gain Phase Analyzer. The four wire measurement 
technique was used to measure resistance R and reactance X. 
The capacitance probe was on one side of the DUT, tank 


water, and the other side was the Harp or a paper clip. 
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Figure 3.14. Capacitance vs. Relative immersion for T-15 
sino precision blocks (B2P108). 
The capacitance can be determined for the measured 


reactance, as a function of frequency, [Ref. 31] by 
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The probe wire was mounted on the Harp and immersed 
approximately 25.0 + 0.5 cm. The theoretical capacitancemmem 
a cylinder of this length is 60.0 pF. The measured 
capacitance was 62.5 + 0.1 pF at 30 kHz using the eee 
clip, and 66.5 + 0.1 pF using the Harp, as shown in Figures 
(3215 and Season 

The static capacitance sensitivities for the various 
measurement techniques are shown in Table (3.2). The 
measured sensitivities are in general good agreement with 


theory. 
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Figure 3.15. R £ X vs. f for Memo wire (T-15) in the wind- 


wave tank using a paper clip (B2P93). 
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TABLE 3.2. STATIC CAPACITANCE SENSITIVITIES FORM 


MEASUREMENT TECHNIQUES. 


Method Sensitivity 
(+0.1 pF/cm) 
T-13 immersed in graduated 
cylinder with tank water at 30 
kHz (B2P23) 
T-15 immersed A graduated 
cylinder with tank water at 30 
kHz DEZ 
Omega TF-101 ] Mercury 
(B2P102) 
T-15 immersed on Harp in wind- 
wave tank using precision 
blocks iba bes 
T-15 immersed on Harp in wind- 
wave tank at 30 kHz, Harp 
ground (BABA) 


T-15 immersed on Harp in wind- 


wave tank at 30 kHz, paper clip 





qgaround (B2P93) 
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2. Conductance Probe 

The idea of a four wire measurement was motivated by 
continued electrical grounding problems in the tank when 
using the capacitance probe and the inability to accurately 
measure the electrical properties of water. The four wire 
technique was described in chapter II-D. It was not obvious 
however, what quantity should be measured. This section will 
describe how we decided to measured conductance and the 
selection process for the best conductance probe. 

The type of measurement was the first item to be 
determined. The static calibration technique discussed in 
chapter III-C-2 was used for these test. Tank water was 
placed in either a 1000 ml Pyrex™ graduated cylinder or a 
four liter desiccator. The BNC cables were connected to a 
HP-4194A Impedance/Gain-Phase Analyzer in the four wire 
measurement mode. A Think Jet printer was connected to the 
analyzer by HBIB bus for data recording. 

A (11.5 x 43.0 cm) (1 x w) Vector perforated printed 
Circuit board had 20, number 43 holes, drilled 0.125 inches 
from one end. The holes were 0.0125 inches apart so that > 
inch, 2-56 stainless steel machine screws could be placed in 
the board. These machine screws became the restraining 
mechanism and electrical connectors for the various samples. 
Two BNC cables were tied to the Vector board. Two Kings, KC- 
79-35, E-Z mini hooks were then attached to the machine 
screws and BNC cables. The Vector board was clamped to a 
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Starret Co., No. 254, master vernier height gage with a one 


inch C-clamp. The experimental setup is shown is Figure 


(3. 1085. 


Think Jet Printer [RSS 
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Graduated Cylinder 


Tank Water 


Figure 3.17. Conductance experimental test setup. 


Resistance R, reactance X, conductance ER and 


Susceptance B, were all measured as a function of wire 


immersion. A 30 kHz, one volt signal was used to drive the 


four 0.0195 inch diameter brass wires that were 0.125 inches 


apart. Brass was used only because it was most accessible 


at the time. Conductance proved to be linear as shown in 


Figure (3.18), which is the type of measurement we want. 


The next test, using the same equipment setup, was 


conducted to determine what kind of metal and sensor 


geometry was best, if there was a frequency dependence, and 


if there was any dependence on drive signal or time. 
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Using a swept frequency in two ranges, 100 Hz to 50 
kHz and 100 Hz to 2 MHz, R, X, G, and B were measured as a 
function of depth for each probe. This provided a check for 
any frequency dependence in the measurement. Typical test 
measurement data is shown in Figures (3.19 & 3.20). The x 


and y axis are both linear in all cases. There is no 


frequency dependence evident. 


® Conductance G o 
e 
80 O Resistance R e 


Conductance (mS) or 
Resistance (Ohm) 





Rele ive immer sio ni cm) 


Figure 3.18. Brass wire G £ R vs. Relative immersion 
(B2 BIT8). 

Brass, gold, platinum, and stainless steel were 
tested in a V shaped or cylindrical geometry. The brass were 
mo ea em x 0.288" x 0.005")(1 x w x t) V shaped strips and 
NENI inch diameter, 10.4 cm long wires. The gold was 0.20 


inch diameter, 8.9 cm long wires. The platinum was (16.5 x 
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0.68 x 0.032 em) (1 4 2 shaped strips. The stainless 


steel was 0.20 inch diameter, 10.2 cm long Wires. 
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Figure 3.19. R € X for a stain esca r eee at a 6.0 cm 


relative immersion vs. Frequency (100-50,000 Hz). 
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Figure 3.20. G € B for a stainless steel wire at 


relative immersion vs. Frequency (100-50,000 


59 


7 
| Da a -- = 
> 


B o MKR 30 040.000 Hz 
OH ms REAL SAS ms 
A Se MS 


EA 
kla 


ao U CM 


The stainless steel wire had the highest sensitivity 
of the metals tested, at 8.4 = The brass wire had a higher 
sensitivity than the brass strip. We decided to use 


stainless steel wire because of its high sensitivity. A plot 


of the measurement data is shown in Figure (3.21). 
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Figure 3.21. Conductance of four metals and two geome comain 
vs. Relative immersion (B2P114-119). 


Figure (3.22) shows that there is no dependence on 
oscillator drive for brass or stainless steel wires, at a 
constant immersion depth of approximately seven centimeters. 

The final test for the conductance probe wire 
selection was to determine if there was a time dependence 
for the stainless steel wire. There is a very small 


mS 


dependence of Oe OR: shown in Figure Na 
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Figure 3.22. Conductance of two metals vs. Oscillator drive 


(B2P148) . 
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Figure 3.23. Conductance of stainless steel wire vs. Time 


FEZ 
61 


A five wire configuration of stainless stec 
was tested and selected for use in the conductance method. 
This configuration was selected for the probe because two 
separate two wlre measurements could be made simultaneously 
at approximately the same location. It was anticipated that 
the two outputs might be used in a gradient measurement. 

The probe has five, 0.020 inch diameter, five inch 
long stainless steel wires that are held at the top by a 
(1.5 x 1 inch) Vector board, and held at the bottom 
(1.5 x 3/8 inch) double sided copper printed circuit Boda 
(PCB). Five 0.020 inch diameter holes were drilled 1/4 inch 
apart, in a line, on both the PCB and the Vector board. 

A Dremel® Moto-tool®, model 395, was used to cut 
electrical pads for each hole on both sides of the PCB. The 
Stainless steel wires were inserted in the holes of each 
board. The wires were then soldered on each side of the PCB. 
A 1/4 inch diameter semi-circular loop was formed at the top 
end of the wire, and one end was soldered into a Vector 
board hole along with a 20 inch long, 22 gage, wire used to 
make electrical connection to the associated electronics. 
The loop was used to adjust the tension of each individual 
stainless steel wire. 

A 2.5 inch diameter "J" was formed on one ends 
20 inch long, 1/8 inch diameter, brass rod. The short SMS 


Of the J was inserted into a Wane inch hole in aa 
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electrically isolated pad on the PCB and soldered on both 
sides. 

A two inch right sided triangular sheet of brass was 
soldered to the other end of the brass rod. A two inch long, 
8-32, brass bolt was soldered to the opposite side of the 
triangular sheet. The screw was inserted into a 1/4 thick 
piece of Plexiglas which supported the Vector board. The 
Plexiglas gave the top of the probe some rigidity. A 8-32 
brass nut was used on the bolt to adjust the tension on all 
five wires. 

The brass rod was bolted to a 3/4 inch diameter PVC 
pipe. The 20 inch wires were fed through the PVC. The PVC 
could then be clamped and the conductance electronics could 
be mounted to the other end of the PVC. The five wire probe 


ME shown in Figure (3.24). 


——20 Gage Wire 
AS 








Plexiglas ~. 
Vector Board—| 


Stainless Steel Wire— 
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Etre DA. Five wire probe construction. 
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E. MEASUREMENT SYSTEM ELECTRONICS 

The third experimental component is the associated 
electronics, which will convert whatever the sensor is 
measuring into a usable electric signal. This component must 
also be linear, sensitive, and immune to noise. The 
constraints and design for the electronics, test procedures, 
results, and final electronics selection will be discussed 
in this section: 

1. Capacitance System Electronics 

The capacitance system electronics is a two channel, 
capacitance bridge, op-amp circuit. The purpose of this 
circuit is to provide an electrical output for the Seema 
probe by converting the capacitance changes of the probe 
wire into voltage changes which could be. iinpur 
conventional electronic signal processing instrumentation. 
The two sections of the circuit must combine to give al 
sensitive and linear output which represents the height of 
the wave on the sensor probe. 

The first section of the circuits aroma 
transformer fed capacitance bridge, as shown in Figure 
(3.25), which provides a means to detect small signals. A 
lock-in amplifier oscillator output is fed into the bridge 
via transformer Tl. The transformers are manufactured by 
ADC, model number TF4RX17YY. They have a frequency range 


from 40 f2ºto 100" KER 
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The sensor probe and water ground form one leg of 
the bridge, variable capacitor Cprobe. The other leg has 
three capacitors in parallel, one fixed value Cl, and two 
variable value. There is a rough adjust C5 and a fine adjust 
C4. The variable capacitors are used to electrically null 
the bridge when the sensor is in the tank water and the 
water is in an equilibrium state. The bridge can also be 
nulled when a dummy load capacitor is inserted into the 
circuit instead of the tank water. 

The typical equilibrium capacitance of the sensor 
probe 1S approximately 60 pF. The other leg of the 
capacitance bridge is the sum of the fixed, rough adjust, 
and fine adjust capacitors. The rough adjust capacitor has a 
value of 27 + 13 pF. The fine adjust capacitor has a value 
of 10 + 5 pF. The fixed capacitor in the adjustment leg is 
then 


Cfixed = Ssensor ` Crough 7 Cfine: 
Ban. 3.4 


The approximate value for eee SIZE. 

The second section of the circuit provides output 
impedance matching and signal gain using an op-amp. The OPA- 
111 op-amp is used with a gain adjust capacitor C6. The OPA- 
111 provides a typical open-loop gain of 120 dB and typical 
common-mode rejection of 100 dB. Complete characteristics 
are found in the manufacture's specification sheet. [Ref. 


22 ] 
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The two channels are identical except for the ground 
used after the transformer. The different grounds arc SER 
for circuit isolation and to minimize the channel emi 
talk. The two channel outputs are matched and maximized 
using the gain adjust capacitor C6, in the op-amp feedback 
cinc LE: 

The electronics had one linearity test. Linearity is 
the most important feature required for the dynamic wave 
height measurement system. The principle of superposition 
can be used if the sensor probe and electronic system 
components each have linear responses. The electronic system 
1S measured using "dry" techniques. 

The "dry" calibration of the electronic system isa 
lab bench experiment using mylar capacitors as a substitute 
for probe wire water capacitance. The values of these 
capacitors were all measured at room temperature using a HP- 
4192A Impedance Analyzer. The nine capacitors used for the 
dummy load test were measured and the results of this test 
are shown in Table (3.3). 

A EG&G Princeton Applied Research, model 5210, 
Lock-In Amplifier was used as a signal source and detector. 
A 31.3 kHz, 1.999 volt signal was applied to the input of 
the electronic system. When both channels are being used, 


two lock-in amplifiers were used. 
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Figure 3.25. Capacitance system electronics circuit. 
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ES 


TABLE 3.3. DUMMY LOAD CAPACITANCE VALUES (B2P16). 


Measured Capacitance 
(+ 0.1 pr) 





Different prime numbers are used for the frequencies 
of both lock-in amplifiers. We used 28.3 and 31.3 kHz 
throughout this experiment [Ref. 33]. 

The output of the electronic system was applied to 
the input of the lock-in amplifier. The typical TOCAN 
amplifier setup configuration is shown in Table (3.4). A 
Kikisul, COS 5060, dual channel oscilloscope was used to 
monitor the electronics system input and output signals. 
Figure (3.26) shows a block diagram of the experimental 
setup. 


TABLE 3.4. LOCK-IN AMPLIFIER SETUP CONFIGURATION: 










| Control Button UREA 
Tracking 


L 
Resolution 
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COS-5060A 





Figure 3.26. Dry electronic system experimental setup. 

The first step of the dry calibration is to zero the 
bridge. The variable capacitor C4 and C5 are adjusted to 
zero the output voltage when a 54.0 pF dummy load capacitor 
is placed in each circuit. 

The second step in the dry calibration is to match 
and maximize both channel outputs. Variable capacitor C6 is 
adjusted on both channels when a 10 pF dummy load capacitor 
1s moved back and forth between the two channels. The output 
from a 450 pF range of capacitors is then measured to test 
the linearity of the electronic system. 


Channel one sensitivity was measured at 3.21 + 0.02 


è aa V 
Sl Channel two sensitivity was measured at 3.41 + 0.03 a 


The linearity for each channel is illustrated in Figure 


me). 


GJ 


1400 


y 
x Ch 1 
1200 X Ch2 


Output vottage (mV) 600 


Capacitance (pF) 


Figure 3.27. Capacitance system electronics output voltage 
vs. Dummy load capacitance (B2P17). 
2. Conductance System Electronics 

The five wire probe illustrated in Figure (3.24) is 
connected to two lock-in amplifiers to make a four wire 
water surface height resistance measurement, as shown in 
Figure (3.28). 

Conductance is inversely proportional to the 
magnitude of the complex impedance, or AV, but is 
proportional to the water height on the probe as described 
earlier in Chapter II (C) (3). The conductance system 
electronics is a two channel division circuit. The purpose 
of this circuit is to invert and balance the two lock-in 


amplifier output signals. The two channels are identical 
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because we thought we might perform a gradient measurement 


using the difference outputs. Gradient measurements were not 


completed. 
PWR SUP 
Lock-in Amplifiers 
| a 
HA ch1 Chi 
In Out q In 
Ch2 Ch2 
A B Out Ref 
Conductance 
Electronics 
Chi 
[xn 
Ch2 
123 45 Oscilloscope 
Probe 


Figure 3.28. Five wire probe electronics setup. 

The output voltage of a lock-in amplifier used in 
the differential mode, provides the input to an Analog 
Devices AD 534 Internally Trimmed Precision IC Multiplier, 
as shown in Figure (3.29). The output of the divider is 


determined by [Ref. 34]: 


10 Zo 


Vout = X4 taby 





Fan Es 
where x1 1s the signal input, yq] is the summing input, and 


Z2 is the numerator input. 


T] 


Resistors Rj, through Rs form a voltage divider 
network which determines the overall gain of both channels. 
Resistor R3 1s adjusted to equalize the two channel outputs 
for the same input signal. The circuit is powered by a KEPCO 


Multiple Output Power Supply, MPS 620M. 
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Figure 3.29. Conductance system electronics circuit. 
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F. WIND-WAVE TANK EXPERIMENTS 
1. Introduction 

There were five measurements made in the wave-wind 
tank. They were the wave power spectrum, surface gravity 
wave phase speed, surface gravity wave wavelength, and wave- 
wind tank system free decay Q and impedance measurements. 
This section will detail the apparatus and techniques used 
for each measurement. 

2. Power Spectrum Measurement 

All the wave-wind tank measurements except the tank 
impedance were based on the output of the power spectrum 
measurements as described in earlier sections. A signal was 
detected using either the capacitance or conductive 
measuring technique and then power spectrum measurements 
were made using a HP-3562A Dynamic Signal Analyzer. 
Measurements were made at various positions in the tank, and 
various wave-wind states. The basic states are: 

a) baseline noise, 

b) various wind speeds, 

c) various driven surface gravity wave speeds, and 

d) combinations of wind speeds and driven surface 
gravity wave speeds. 

Data was recorded using a HP-7470A plotter or stored 
on disk using a HP-9133 hard disk. To store data on disk, 
the HP-3562A front panel control sequence in Appendix C was 
used. A HP-54501A, 100 MHZ digitizing oscilloscope was used 


To 


to monitor the different electrical output signals RA 


experimental electronics setup is shown in Figure (3.30). 
HP-S4501A 





HPIB HP- 74702 
Bus LIC | 








neato 


UL UU LL LL Le 


Figure 3.30. Measurement electronics setup. 

Power spectrum measurements were taken with the HP- 
3562A using power spectrum measurement and display modes. 
Tvpical HP-3562A setup configurations are shown in Figure 


Dudo. 
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Figure 3.31. Typical HP-3562A power spectrum setup 
configuration. 
3. Phase Speed and Wavelength Measurements 

The phase speed and wavelength were measured using 
the cross-correlation output of the HP-3562A Dynamic 
Analyzer. Two signals had to be detected to make the phase 
speed measurement as described in chapter II-F using the 
experimental setup shown in Figure (3.30). 

The HP-3562A computed and displayed the cross- 
correlation of the two signals. Typical HP-3562A setup 
configurations are shown in Figure (3.32). The Phase speed 
was then computed using the method of least squares fit on 


T9 


the data output of the cross-correlation measurement. The 


distance between the two probes had to be measured, then the 


wavelength could be computed as described in chapter II-B. 
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Figure 3.32. Typical HP-3562A cross correlation setup 


configuration. 


4. Free Decay Q Measurement 


The free decay Q was calculated using the power 


spectrum measurement output of the HP-3562A over a period of 


ten minutes. 


The ten spectra were measured one minute apart 


and stored on disk for each data set using the experimental 


setup shown in Figure (3.30), as described in chapter II-H. 
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The entire data set was latter plotted on the same 
paper. The data required for any measured frequency could be 
extracted from the data on the disk files, and the Q could 
be calculated. 

5. Tank Electrical Impedance Measurement 

A self contained mini experiment was performed to 
answer the question of why the lock-in amplifier oscillator 
was being loaded down by 5-50% during various capacitance 
measurements. The tank impedance, via the water and 
plumbing, returning to the rack ground was tested. 

The lock-in amplifier output voltage was measured on 
an oscilloscope with and without the tank impedance in the 
Circuit using the equipment setup shown in Figure (3.33). 
The tank impedance is then calculated using the measured 


voltage at any given frequency, by 


Bam. 3.6 
where Zi] 1s the lock-in amplifier output impedance. The 
lock-in amplifier has a 600 ohm output impedance which was 


verified with a 100 ohm resistor. 
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Figure 3.33. Tank electrical impedance experimental setup. 
Voltage output for the closed circuit as a funda 
of frequency is plotted and shown in Figure (3.34). The open 
circuit voltage was constant at 5.7 Vop: The measured tank 


impedanceyat 20 kHz is Z(20 KHZ) = MSN: 
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Figure 3.34. Lock-in amplifier output voltage vs. Frecuemas 


for tank BA 


78 


IV. RESULTS 


A. STATIC CALIBRATIONS 

The static calibrations were certainly the easiest 
measurement of the experiment. They were easy because there 
was no complicated equipment, the technique described in 
chapter III-C-2 is very easy to follow, and we used the 
static calibration technique so many times. 

1. Capacitance Static Calibrations 

All of the capacitance probe static calibrations 
were in good agreement with the theoretical capacitance 
Peioitivity of 2.4 + 0.1 pF/cm as shown in Table (3.2). 

A theoretical capacitance system static voltage 
sensitivity could be determined by combining either the 
theoretical capacitance sensitivity or the measured Teflon 
wire capacitance sensitivity with the measured electrical 
sensitivities of each channel of the system electronics. The 
theoretical voltage sensitivities for each channel are shown 
in Table (4.1) using the results from the previous chapter. 

Capacitance system static voltage sensitivities were 
measured using the static calibration technique with water 
contained in the Pyrex graduated cylinder and the 
electronics shown in Figure (3.26). The lock-in amplifier 
oscillator was set at 30 kHz for these measurements. A plot 


of the measurement data is shown in Figure (4.1). 
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TABLE 4.1. THEORETICAL CAPACITANCE SYSTEM STATIC ACTE 


ITEM SENSITIVITIES - 
Ee 
ac E 
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Figure 4.1. Capacitance system electrical output voltage vs. 
Relative immersion in the graduated cylinder (B2P24). 

Capacitance system static calibrations were also 
conducted using the Cone Head apparatus described in 
Appendix B. Water was added to the Cone Head before each 
measurement during this test, using Omega thermocouple wire 
probe T-13. A plot of the measurement data is shown in 
Figure (4.2). 

The capacitance system static voltage sensitivities 
for channel one measured in the graduated cylinder are in 
good agreement with theory. It is unknown why the 
capacitance system voltage sensitivities for channel two 
measured in the graduated cylinder were not in better 


agreement with theory. 
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Figure 4.2. Capacitance system electrical output volta ON 
Relative immersion in the Cone Head apparatus (B3P2). 

The system voltage sensitivities measured using the 
Cone Head apparatus are not in good agreement with theory. 
This is probably because of the apparatus and procedure 
itself. We could accurately measure the water that was being 
put into the Cone Head system. We could not accurately 
measure the amount of water that was actually getting into 
the changing water column. Water would stick on the sides of 
the graduated cylinder and the wire. 

The theoretical and measured sensitivities for all 
the test, and the percent differences are shown in Table 


(4.2). 
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2. Conductance Static Calibrations 







The properties of ordinary water are so complex that 
a theoretical conductance sensitivity 1S impractical. The 
important thing for these calibrations was not how do they 
compare with theory, but are the calibrations linear, 
repeatable, and measurable. 

Several wire conductance static calibrations were 
made using the static calibration technique discussed in 
chapter III-C-2. One series of static measurements was made 


using the HP-4194A Impedance\Gain-Phase Analyzer with the 
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wires in either a ten gallon Pyrex beaker or the wind-wave 
tank. These measurements were made to compare the results 
from data taken in the two water samples. 

The conductance sensitivity of the measurements 
taken in the Pyrex is 6.4 + 0.1 mS/cm. The conductance 
sensitivity of the measurements taken in the wind-wave tank 


is 6.3 + 0.1 mS/cm. The data from both measurements RS 


linear, as shown in Figure (4.3). 
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Figure 4.3. Conductance vs. Relative immersion in two water 
systems (B3P31). 

Another series of static wire measurements were made 

using the same equipment, in the wind-wave tank. These 

measurements were made to see if two four-wire measurement 


combinations were possible using the five wire probe, and 
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how close the sensitivities were for the combinations. A 
gradient measurement technique was discussed, but not 
pursued. Table (4.3) shows the possible combinations. The 
wires that are used in each combination are shown filled in. 
MABLE 403. FOUR WIRE COMBINATIONS WITHIN THE FIVE WIRE PROBE 


(B2 P3 9i 


Four Wire Combination Four Wire Combinations 


Number 





We used deviations in the conductance measurements 
to ensure that all the individual wires were straight as 
possible. We did use the tensioning loops to adjust some of 
the wires. The final measured conductance sensitivity and 
deviation results are shown in Figures (4.4 and 4.5). 

The conductance sensitivities, shown in Table (4.4), 
of the various four-wire combinations are different because 
of the two different probe geometries, four in a row or a 
space inserted. These geometries put the voltage probes in 
different potential fields between the current probes. The 


sul and SU3 configurations have lower sensitivities than the 


SU2 and SU4 configurations because the AV is the lowest in 
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this symetric configuration. The highest sensitivity should 
be obtained when the center probe is not used. 

The conductance sensitivity measurements were 
linear, with high correlation Cocftfireremes and the 
deviations were consistent for all the combinations with a 
maximum error of 0.4%. The causes of the anomalies of the 
last data points for SUl and SU4 are unknown. We choose to 
use wire combinations SU2 and SU4 in all latter test because 


they had the highest conductance sensitivities, as expected. 
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Figure 4.4. Conductance vs. Relative immersion for various 


four wire combinations (BOPA 
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Figure 4.5. Deviation vs. Relative immersion for various 
four wire combinations (B2P42). 
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Conductance system static calibrations were made to 
determine the voltage sensitivities which would become 


inputs for the engineering units of the HP-3562A Dynamic 


87 


Analyzer power spectrum measurements. The lock-in amplifiers 
provided electrical inputs to the conductance electronics. 
The conductance electronics electrical output was then 
measured using two hand held multimeters, one for each 
channel. 

The voltage sensitivity for channel one is 0.309 + 
0.005 V/cm, and for channel two is 0.381 + 0.005 VCCI RS 
data from the measurements was linear, as shown in Figure 


(4.6). 
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Figure 4.6. Output voltage vs. Relative immersion for 


conductance system electronics (B3P59). 
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B. DYNAMIC CALIBRATIONS 
1. Capacitance Dynamic Calibrations 

Several capacitance dynamic calibration apparatus 
and techniques were tried, as described in Appendix B. The 
apparatus was cumbersome, and the results of the techniques 
were not repeatable. The repeatability problem was caused by 
the low tank impedance, the results of this test were 
described in the previous chapter. 

It was also determined that all the capacitance 
dynamic calibration techniques were not conducted in the 
preferred motion either. These techniques moved the wire up 
and down in the water. This does not represent a traveling 
gravity wave motion. 

The low, approximately 20 ohm, tank impedance made a 
two wire capacitance system unusable when combined with the 
unknown and unstable water properties. Unfortunately, the 
tank impedance test was not conducted until late in the 
Capacitance data collection efforts. 

2. Conductance Dynamic Calibrations 

Conductance dynamic calibrations were conducted to 
determine repeatability, if there was a preferred probe 
orientation, and to determine the frequency response. The 
calibrations were conducted in the small wave guide to cut 
off low frequency parametric waves, using the techniques 


described in chapter III-C-3. 
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Four consecutive measurements were taken in the same 
location, with three different probe orientations. Thee 
of wires in the probe were either normal or parallel to the 
paddle. The first measurement of the sequence was always 
parallel to the paddle. Two measurements were always made 
with the wires parallel to the paddle. 

The two parallel measurements were the test of 
repeatability. The other measurements were the test for best 
probe orientation. The probe frequency response could be 
determined at all orientations. 

A two Hertz, one to four Hertz, dynamic power 
spectrum calibration was conducted in the small wave guide. 
Waves below two Hertz were not sinusoidal on the 
oscilloscope, which set the lower frequency limit. The Flat 
Top window was used for these measurements because this 
window has the lowest equipment, HP-3562A Dynamic Analyzer, 
amplitude accuracy uncertainty of up to + 0.1 dB [Ref BENIN 

The engineering units used in the HP-3562A internal 
power spectrum calculations were determined from the results 
of the static calibrations discussed in section A(2). 

The parallel probe orientation consistently had the 
highest sensitivity at any given frequency for both 
channels, as shown in Figure (4.7 and 4.8). The data pa 
at four Hertz is above the parametric wave cutoff for the 
small guide. The last daa points show a definite 
interaction between the gravity and parametric waves for all 
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probe orientations. The parallel probe orientation was also 
very repeatable, except for the anomaly at the 2.5 Hertz 
knee. 

The HP-3562A Dynamic Analyzer measurements were 
verified by measuring wave peak-to-peak amplitudes using a 
HP-54501A Digitizing Oscilloscope. The amplitudes ranged 
mom 330 mV at two Hertz to 187 mV at four Hertz (B3P68-90). 

The very disappointing result was that the probe 
frequency response was not flat, not even for an octave. We 
have assumed that the probe response would continue to drop 
with frequency, but could not test the higher frequencies 
because of the parametric waves caused by the paddle. The 
current probe construction would not allow for a smaller 
wave guide, which would enable us to test at higher 


frequencies. 
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Figure 4.7. Sensitivity vs. Frequency for channel one 


(BSPO2=67). 
ol 


O 


xe 
ex 


Sensitvity (6B) 








e il 
O Per 
x Il 


x Per(180) 


Frequency (Hz) 


Figure 4.8. Sensitivity vs. Frequency for channel two 
(B3P64-67). 
C. TANK MEASUREMENTS 

Measurements were taken in the wind-wave tank through- 
out the experiment. The sensor may or may not have been 
statically calibrated prior to any particular measurement. 
No tank measurements were made after the conductance system 
dynamic calibrations were made. 

Results of early phase speed and power spectrum 
Capacitance system measurements were presented at the 119th 
meeting of the Acoustical Society of America, at Penn State 
University [Ref. 36]. 

The cross correlation plot in Figure [4.9] is a typical 


capacitance system measurement phase speed plot. The data 
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were analyzed using the least square fit technique. This 
1.35 + 0.01 Hertz, driven gravity wave has a measured 
frequency of 1.36 + 0.01 Hertz. The measured delay time is 
15.6 + 0.1 msec. The probes had a measured separation of 
0025 + 0.01 m. 

The theoretical phase speed for this driven gravity 
REC using equation [2.3], is 1.16 + 0.01 m/sec. The 
measured phase speed, using equation [2.25], is 1.63 + 0.01 
m/sec. This is a 40.7% difference. It is possible that 
electrical cross-talk between the two channels in the 
electronics box may have decreased the delay time measured. 

The theoretical wavelength for this same wave, using 
equation [2.26], is 0.86 + 0.01 m. The measured wavelength 
RS 1.20 + 0.01 m. This is a 39.9% difference. 

The capacitance system equipment and techniques were 
both adjusted, and the measured phase speed and wavelength 
were consistently higher than theory. We believe the results 
are high due to an erroneous delay time in the channel two 
output caused by electrical cross-talk in the electronics 
box, but this was not confirmed. 

Base line noise and wind power spectrum measurement 
plots, using a statically calibrated capacitance system, are 
shown in Figure (4.10). Similar data obtained with the 
conductance system is shown in Figure (4.11). The four upper 
lines in Figure (4.10) are wind power spectrums, at four 
different wind speeds. Spectra for both the capacitance and 
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conductance systems peak at approximately two Hertz and 
decay at approximately 50 to 70 dB per decade, or with power 
decaying at 07% to 077, but the slope of the decay malh not 
be trustworthy because the transducer was not dynamically 
calibrated and its frequency response is unknown. However, 
there is excellent agreement in the wind power spectrum 


measurements for the two transduction systems. 
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Figure 4.9. Capacitance system cross correlation picas 
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Figure 4.10. Capacitance system base line noise and wind 


power spectra plot. 
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Power spectra of (a) Base line noise, (b) driven gravity 
wave, and (c) driven gravity wave with wind, using a 
statically calibrated capacitance system, are shown in 
Figure (4.12). The driven gravity wave spectrum shows the 
harmonics of the driving frequency, at least 18 dB below the 
fundamental. There is very little energy above the noise 
level at other frequencies. 

The driven gravity wave and wind spectrum is 
approximately 20 dB above the noise and the driven gravity 
wave spectrum. The wind energy completely dominates the 
harmonics of the driven gravity wave with this particular 
frequency resolution. The wind decay is unaffected by the 
addition of the driven gravity wave. The decay of the power 
is approximately as we. 

The power spectra plots of the driven gravity wave and 
driven gravity wave with wind were expanded to 0.5 Hz about 
the fundamental driven frequency, shown in Figure (4.13). 
The energy peak is downshifted 0.062 Hz when wind is added. 
A second, lower frequency energy peak 1s also observed. This 
peak 1S approximately 30 dB above the noise level. This peak 
may be caused by the non-linear effects discussed in Dr. 
Larraza's theory. Determining the cause(s) of the additional 
energy peak and the fundamental frequency downshifting were 


Mert for follow-on theses. 
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Figure 4.12. Capacitance system base line noise, driven 
gravity wave, and driven gravity wave with wind power 


spectra plot 
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Figure 4.13. Expanded plot of driven gravity wave and driven 


gravity wave with wind. 
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Plots of the tank quality factor measured with the 
Capacitance system, are shown in Figure (4.14) and similar 
data obtained with the conductance system is shown in Figure 
(4.15). The tank Q measurement with the conductance system 
was almost double the capacitance system measurement 's value 
at the low frequencies, and half the value at the higher 
frequencies. It is unknown why and how there is a difference 
in tank Q for the two measurement techniques. The measured 
Q at different frequencies is shown in Figure (4.16). 

The driven gravity wave spectrum was not measured 
using the conductance system. All efforts were directed at 
getting a dynamic calibration for the conductance system 
instead. We were not sure of the best method to make phase 
speed and wavelength measurements with the conductance probe 


at that time. 
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Figure 4.14. Tank Q capacitance system measurement plot. 
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Figure 4.15. Tank Q 
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conductance system measurement plot. 
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Figure 4.16. Tank quality factor vs. Frequency measured by 


Capacitance and conductance systems (B2P137-142 and B3P10). 
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V. CONCLUSIONS 


Capacitance and conductance measurements using two and 
four wire techniques were developed, tested, and calibrated 
in this thesis. The four wire conductance system technique 
was determined to be the better of the two tested, and is 
currently being used in a follow on thesis. The two wire 
capacitance system repeatability problems caused by 
grounding and the complex properties of "simple" water are 
not present using the four wire technique. The two systems 
did give quite similar results for wind-wave spectra. 

We feel the four-wire technique 1S superior to the two- 
wire technique as long as the potential problems in 
conducting a water surface wave height measurement, dirty 
probes and unknown variable water chemistry, only exist ona 
microscopic scale. This thesis did not prove or set limits 
to this assumption. A study of these problems is suitable 
and desirable as a follow on thesis. 

Static calibrations were performed many times on all 
systems. The results of the various tests were consistent, 
linear, and generally compared well with theory. The static 
calibration techniques are simple and should be completed 
prior to other measurements in future experiments. 

The results of the conductance system dynamic 
calibration were disappointing and puzzling because segs 


frequency response was not flat. It is difficult CONS 
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what the frequency response is in the entire frequency range 
of interest, 0-10 Hz, due to the limited results we were 
able to obtain. However, the two systems gave quite silimar 
results for the wind-wave spectra. 

Driven gravity wave frequency downshifting and wind 
energy dumping was observed in the combined gravity wave and 
wind-wave measurements. A follow-on thesis will explore 
enis. 

The importance of the dynamic calibration is unknown for 
future nonlinear experiments and should be reviewed to 
insure that important data is available and interpretations 
are accurate. 

The most important thing I learned from this research 
are some basic principles of experimental science. First, 
you must have some idea concerning what and why you are 
investigating nature. Second, you must have a predetermined 
purpose and procedure for the investigation. Third, data 
must be recorded thoroughly, meticulously, and in a timely 
manner. Fourth, you must evaluate the data carefully to 
determine the outcome of your experiment. Finally, you need 


GOOD LUCK. 
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APPENDIX A 


I. HP-15C CALCULATOR LEAST SQUARES FIT STANDARD DEVIATIONS 

PROGRAM CODE 

The standard deviations in all the least squares fit 
parameters were discussed in chapter II E. This is program 
code for the HP-15C hand held calculator, which can Dewees 
to solve equations (2.16 through 2.19) for data entered into 
the calculator using the summation key (B1P66/67). The 
calculator stores the different summation values in the 
following register as shown in Table (A.1). 


TABLE A.1. HP-15C REGISTER VALUES AFTER USING SUMMATION a 


N 











| a | summation of xt 
|g omaan of ga AH 


The standard deviations in y, slope, and intercept are 





displayed at program lines 040, 054, and 067 respectfully. 


The program code is shown below. 
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APPENDIX B 


I. CAPACITANCE SYSTEM DYNAMIC CALIBRATIONS 

This section will briefly describe the equipment and 
techniques used for the capacitance system dynamic 
calibrations. Two dynamic calibrations of the capacitance 
system were explored before the capacitance system was 
determined to be unusable. Both of the techniques moved the 
water or the sensor probe, up and down, which we decided was 
not the preferred motion for the dynamic calibration. 

A. Cone Head 

A 16 ounce plastic funnel, or cone, was used to 
connect a 110 ml glass burette to a USRD ‘Type J-11 
Transducer, hence the Cone Head name for the apparatus. The 
J-11 provided a means to convert an electrical signal to 
fluid motion and the funnel acted as a fluid amplifier with 
a gain of approximately 130. All capacitance system 
measurements were actually made in the burette. 

The top of the probe wire was held in place at 
the top of the burette with a small cork stopper. Fishing 
line was tied to the bottom of the capacitance probe. The 
fishing line was passed around a large paper clip lodged in 
the funnel, pulled out of the top of the burette to tension 
the probe wire and held in place at the top of the burette 
with the same small cork stopper. 


EOS 


All joints in the Cone Head were sealed for 
water tight integrity with Polymeric Systems, Inc., PSI 601 
Clear, 100% Silicone, RTV Adhesive Sealant. The Cone Head 


apparatus is shown in Figure (B.1). 


——— 110 ml Burette 


<— 16h oz Funnel 


Jl TEN Sane I 


/ Pressure Valve 


Figure B.1. Cone Head apparatus. 

The J-11 transducer motion was calibrated using 
an LVDT. The J-11 compensation bladder pressure was measured 
using a PX-126, Omega Engineering, Inc. differential 
pressure transducer. This provided a non-intrusive 
measurement of the J-11 when the Cone Head was assembled. 

The electrical outputs of the LVDT and pressure 
transducer were compared when a low frequency square wave 
was applied to the J-11. The Cone Head was assembled with 
water in the apparatus and the electrical output of the 
pressure transducer was then compared with the capacitance 


system output. 
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This calibration scheme was not refined due to 
the fact that only one sensor probe could be fit into the 
Cone Head at a time and the poor results in the Cone Head 
Seatcic Calibrations. 

B. Motor 

This technique was briefly tested, with no 
documented results. The sensor probe was moved up and down 
in the wind-wave tank, while the probe was on the Harp. 

A 1/4 horsepower AC motor was bolted to the top 
U channel on the Harp. A 2.5 inch diameter pulley was 
mounted eccentrically to the motor shaft. A one inch 
diameter pulley was mounted directly above the probe. One 
end of an approximately 24 inch long, ten pound, fishing 
line was tied to the top of the probe, placed in the grooved 
flywheel, and the other end of the fishing line was tied to 
a standoff post on the 2.5 inch pulley. The other end of the 
probe was connected to the lower U channel with two number 
one rubber bands. 

The AC motor was plugged into a variable 
transformer which would adjust the mechanical output 
frequency of the large flywheel and provide the forcing 


Mmomction for the probe. 
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APPENDIX C 


I. HP-3562A DYNAMIC SIGNAL ANALYZER FRONT PANEL CONTROL 
SEQUENCE TO STORE DATA ON DISK 
The front panel hard disk accessing sequence for the HP- 
3562A Dynamic Signal Analyzer using the HP-9133 Disk Drive 
1s shown in Table (C.1). Complete operating instructions can 
be found in Reference 35. 


TABLE C.1. HP-3562A HARD DISK ACCESSING TS EDEN 


AN HP-IP FCTN Select Address 


a SELECT ADDRESS Disc Address = O 
Disc Unic = 


La [orse save/Recall File 
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